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Abstract. Non-Markovianity has recently attracted large interest due to significant advances in its 
characterization and its exploitation for quantum information processing. However, up to now, only non- 
Markovian regimes featuring environment to system hackflow of information (strong non-Markovianity) 
have been experimentally simulated. In this work, using an all-optical setup we simulate and observe the so- 
called weak non-Markovian dynamics. Through full process tomography, we experimentally demonstrate 
that the dynamics of a qubit can be non-Markovian despite an always increasing correlation between 
the system and its environment which, in our case, denotes no information backflow. We also show the 
transition from the weak to the strong regime by changing a single parameter in the environmental state, 
leading us to a better understanding of the fundamental features of non-Markovianity. 


Introduction 

The development of quantum technologies for information processing, communication and high 
resolution metrology among other applications has renewed the interest in a better understanding of 
the dynamics of open quantum systems. The most typical description of an open system evolution is 
that of a Markovian dynamics caused by the memoryless interaction of a given quantum system with its 
environment fT]. On the other hand, strong system-environment interaction, environment correlations 
or initial system-environment correlations may cause memory effects rendering the dynamics non- 
Markovian. Recently, non-Markovian dynamics has become a very trendy topic mainly due to the 
development of new experimental techniques for controlling and manipulating solid state and many 
body systems HJ |3l 13 |5l |6l and also due to its possible applications in information protection and 
processing EiiiiiiMiiiiini- 

Commonly, the evolution of a system is defined as Markovian if the corresponding quantum map 
is divisible in other completely positive maps (from now on CP maps), i.e. A, 2,0 = A, 2 ,(iAfj 0 for all 
^2 > T > 0 |[T3l IT3l . For all the maps that do not satisfy this equality, the corresponding evolution is 
non-Markovian. The conditions for a strict Markovian dynamics are usually very hard to achieve and 
most experiments will present some degree of non-Markovianity. However, it is not always that the 
non-Markovian characteristics of a dynamics can be easily observed. Some non-Markovian processes 
can be identified in ferms of a measurable quanfify such as an increase in fhe disfinguishabilify of 
differenf quanfum sfafes lITSll or of fhe enfanglemenf befween fhe evolving sysfem and an ancilla ifl^ . 
These processes presenf a strong degree of non-Markovianity that has also been called essential non- 
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Markovianity iflTl as opposed to a weak non-Markovianity that requires full process tomography and, 
therefore, is much more difficult to detect. 

It is important to understand the effects caused by such reservoirs in quantum computational systems, 
such as qubits, and to observe these effects in controlled experimental setups in order to pinpoint the 
essential mechanisms behind the different time evolutions generated by them. Hence, very recently, non- 
Markovianity has been investigated in different setups and contexts such as the control of the initial states 
of the environment |[T8l[T9ll . of its interaction with the system EOll^ or combinations of them |[22l . as 
well as the observation of non-Markovian effects in simulated many-body physics ll2^ or in the recovery 
of quantum correlations |[2^l25l . All these experiments are restricted to detecting and/or exploring the 
strong non-Markovianity. 

In this work we carry out an experimental characterization of the transition between weak and 
strong non-Markovianity. More specifically, the term weak non-Markovian will be used for dynamics 
represented by maps that are divisible in positive maps, but not in completely-positive maps. On the 
other hand, strong non-Markovian will be used to maps that are not even divisible in positive maps. 
In particular we adopted full process tomography to observe the weak non-Markovianity dynamics of 
a qubit subjected to the interaction with a correlated environment. Furthermore, through the careful 
preparation of the environment state, the transition is induced by changing a single experimental 
parameter. 

Theoretical model 

We consider a qubit that interacts with an environment from which it is initially decoupled. The 
interaction consists of consecutive collisions each of which can produce three different effects on the 
system: either nothing happens, in which case the identity is applied to or the system is rotated by 7t/2 
around the X or Z axis, undergoing a or Xg = flip. The map that describes the evolution by 

one collision is Aio(-) = Po{') + + {po-,Px,Pz} probabilities that add to one. 

This is a special case of a random unitary qubit evolution. It is also a unital map. 

All the effects we are interested in can be observed with only two collisions described, in our case, 
by the general map A 2 o(-) = 'LmnPmnOnOni{-)OmOn where 0 < < 1, 'LmnPmn = 1 and Oq = 1, 

Ox = and 0^ = 'Ijs- Note that if the collisions are fully independent (hence non-correlated) this 
model is Markovian by construction and can be easily generalized for any number of collisions. The 
state of the system after n collisions is obtained by the concatenation of single collision CP maps: 
Pi(n) = A„oPi(0) = (Aio)"Pi(0). In this case, the joint probabilities of two consecutive flips need to 
respect relations such as = pipj where i,j = {x,z}- 

The dynamics becomes more interesting if the collisions are correlated, i.e. when pij ^ PtPj- In 
particular, it is shown in Ref. [26] that for any correlation factor Q = larger than zero, the 

two-collision map A 20 represents a non-Markovian evolution, i.e. A 20 / A 21 A 10 or, equivalently, Pi(2) 
cannot be obtained by applying a CP map on Pi(l). Naturally, larger values of Q produce a more intense 
non-Markovian effect. There is, however, a transition in the type of non-Markovianity that depends 
on the probabilities of the flips. If correlated flips are very likely, i.e. pu (i = {x,z}) is of the same 
order of poi, p,o and poo> then the non-Markovianity is strong (sometimes referred to as “essential” in 
the literature ifTTl l which means it can be witnessed by quantities such as the trace distance between 
different states of the system or the entanglement between the system and an ancilla. This is related to 
the fact that the reconstructed map A 21 = A 2 oAjq^ is not even positive, let alone CP, i.e. it does not map 
the Bloch ball onto a set contained in it. The extreme scenario has pxx = = 1/2 in which case after 

the first collision will be given by Pi(l) = 5 Xpi( 0 )X -|- 2 Zpi( 0 )Z and after the second collision it goes 
back to Pi(0) (since X^ = Z? = 1). 

As the flip probabilities decrease this effect becomes smaller and at some point the reconstructed 
map Ai 2 becomes positive, albeit still non-CP. For random unitary maps, as it is the case here and as 
it is shown in Ref. [27], a map is divisible in positive maps if and only if the von Neumann entropy 
and the trace distance present a monotonic decay in time, establishing a strict relation between backflow 
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Figure 1. (a) Conceptual scheme of the experiment. A maximally entangled state between the qubit 
of interest (s) and an ancillary state {a) is initially prepared at time to- The interaction between qubit 
s and the environment is simulated by a sequence of two channels, performing each a mixture of 1, 
X and Z operations. The state after the interaction channel CH\ {CH 2 ) is measured at ti { 12 ). In the 
actual experiment the two channels are simulated each one by two liquid crystal modulators. Qubit 
a will not suffer any change, since it is isolated from the environment, (b) Sequence of probabilities 
corresponding to the action of the two channels, poo = (1 ~ 2s)^, pojc = POz = PzO = PxO = (1 ~ 2s)s, 
and Pjcjc = = 2s^. X; (Z,) is the X (Z) operation occurring in the /th-collision. Here we adopted the 

specific value of s = 0.2. 


of information in terms of these quantities and strong non-Markovianity. In the case of our model, 
this backflow of information can also be identified observing the entanglement between system and an 
ancilla. In this case, the previously discussed witnesses fail and only full process tomography of the 
dynamics at each step can detect the non-Markovianity character of the evolution. This situation has 
been defined as “weak” non-Markovianity in Ref. [17]. The goal of this work is to observe this weak 
non-Markovianity as well as the transition of non-Markovian regimes as a function of the correlated flips 
probabilities. Finally, note that orthogonal flips (X or Z) are chosen to maximize the effect but, in fact, 
any pair of non-commuting flips will also produce a non-Markovian map for 2 > 0 |[26l . To implement a 
Kraus operator for the two collisions we exploit a mixture in time of the different Pauli operators (since 
the time emission of the photon is random). In our scheme we simulate non-Markovian dynamics with 
a classical apparatus; for more details refer to Ref. [25]. The conceptual scheme of the experiment we 
are going to show is depicted in Fig. [T^). Here CH\ and CH 2 are the interaction channels acting on the 
system. Fig. [^) shows the probabilities associated to the sequence of operations performed by CH\ and 
CH2. 


Experimental setup 

In the experiment, the system s is the polarization state of an initial maximally entangled photon 
pair, |\|/ai) = {\HV) + e'^\VH))/x/l generated by a PPKTP ultrabright source of polarization entangled 
photons (281, where \H) (|F)) represents the horizontal (vertical) polarization (see Fig. for details) 
(281 . Note that |/)a ® \j)s will be simply represented as |/ j). The environment is simulated by a sequence 
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Figure 2. Detailed scheme of the experiment. Twin photons are created by a polarization entanglement 
source. One photon (system s) is sent through a correlated liquid crystal environment, while the other 
(ancilla a) is let to go free. Then, the bipartite state is measured by complete state tomography at times 
to, t\ and t 2 - Liquid crystals (LCs) (two for CHi and two for CH 2 ) act as phase retarders, with the relative 
phase between the ordinary and extraordinary radiation components depending on the applied voltage. 


of four voltage controlled liquid crystal cells (LC) lying on the path of photon 5 . By a suitable control 
of the applied voltage on each of them, the four LCs were set to operate either as the identity or as 
half-wave plates. In particular the first and third LCs were oriented with the slow axis along the vertical 
direction, thus acting either as 1 or Z, while the slow axis of the second and fourth LCs were oriented 
along the diagonal direction (45°), thus acting either as 1 or X. According to the collision model, 
the first environment (giving Pai(l) as output result) derives from the actuation of the first two LCs 
(1 and 2) only, leaving the other LCs (3 and 4) in the identity regime. Finally, the second collision 
environment (giving Pai(2) as output result) corresponds to the actuation of the four LCs in the designed 
way. The parameter s gives the probability that either X or Z occurred in the first collision {px = E 
and = e) and it is proportional to the time of application of the voltage to the liquid crystal. In the 
experimental setup it is defined by fhe ratio of the width of an applied voltage pulse to the width of a 
measurement cycle. Since the photons in our source are generated at random, the action of the liquid 
crystal cells driven by the voltage pulses will in fact simulate random collisions. By controlling the 
time duration of the applied voltage on each LC intercepting photon s, it was possible to choose the 
probability corresponding respectively to the 1, X and Z operations. This is given by the parameter 
E in the following way: poo = (1 - 2e)2, pQx = poz = PzO = PxO = (1 - 2e)e, and pxx = Pzz = 2e^, as 
shown in Fig.[^). Note also that pxz = Pzx = 0 reflecting the fact that only perfectly correlated rotations 
(2 = 1) are implemented. In this case, the theory predicts an always non-Markovian dynamics for any 
value of E > 0. In order to verify the dynamical behavior, the measurements were performed after the 
first collision {t\) (in that case LCs 3 and 4 were acting as the identity and only LCs 1 and 2 were varied) 
and the second collision (I 2 ) (with all four LCs varied). As described below, in the case of strong non- 
Markovianity, we need just to measure an entanglement witness. However, for weak non-Markovianity, 
quantum state tomography should be realized. 
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Figure 3. (a) The negative eigenvalue of i?/ as a function of s. The inset is the same curve for small values 
of s, £ < 0.1. For qualitative reasons, the weak non-Markovian regime, where the intermediate map is 
positive, but not completely positive, is represented by the blue region and the strong non-Markovian 
regime, where the intermediate map is not even positive, by the red region. The experimental error 
bars are estimated from propagation of the Poissonian statistics of photon coincidence countings in the 
tomographic reconstruction of the process matrix, (b) The difference between the concurrences of system 
and ancilla after two collisions C(2) and one collision C(l) versus £. The inset shows C(2) and C(l) versus 
£ plotted separately. For qualitative reasons, the weak non-Markovian regime is represented by the blue 
region and the strong non-Markovian regime by the red region. The transition from one region to the 
other is represented by the dashed line, which can vary its position depending on the imperfections in 
the experiment. The experimental error bars are estimated as explained before. The uncertainties, about 
l%-3% of the concurrence values, are within the size of the symbols. 


In the experiment the open system dynamics is obtained by temporally mixing the three possible 
settings of CH\ (1, X, or Z), giving Pai(l), and by temporally mixing the seven possible settings of 
the action of CHi and CH 2 (namely IX, XI, IZ, Zl, XX, ZZ, 11), giving Pai(2). Note that only 
correlated rotations are implemented; this is done in order to maximize the non-Markovian effect as it is 
better explained in Ref. [26]. The map Aio (A 20 ) is obtained from the full tomographic reconstruction 
of Pai(l) (Pai(2)) and the intermediate map A 21 that tells us about the character of the dynamics is 
calculated from A 21 = A 2 oAjq'. 

Non-Markovian analysis 

The density matrix of a qubit state can be represented by p = (1 -|-r • a) /2, where f is the Bloch vector 
(r,- = Tr(pa,)). The action of a map A on p can be described, in general, by A : r 1 —)• r' = Mr+7, where M 
is a matrix responsible for changing the norm and rotating the Bloch vector while 7 = {tx,ty,t^) shifts its 
origin. For unital maps (7 = 0), which is the case studied here, one can define a 4 x 4 Hermitian matrix 

= {\ + <8) CJy) /2 SO that the map A is completely positive iff > 0 1^1311 . 

In Fig. [^) we plot the minimum eigenvalue X of J-C for the map A 21 obtained from the 
experimental tomographic reconstruction of the system state after one and two collisions (Pi(l) and 
Pi(2) respectively). The fact that X is always negative for £ > 0 necessarily implies that A 21 is non- 
CR As a consequence, the dynamics of the system is non-Markovian for any value in this range. The 
question remains whether the measured non-Markovianity is weak or strong. A linear map is positive 
iff the corresponding dynamical matrix H is block-positive |[30ll . For the type of map implemented in 
our experiment, a simple calculation (see Methods) shows that the condition that guarantees the strong 
non-Markovian regime, which means that H is not block-positive, necessarily implies the recovery 






























of entanglement between the system affected by the environment and an ancilla used to monitor the 
dynamics. Therefore, in order to search for an eventual transition from weak to strong non-Markovian 
regime in our case, we have also measured the variation of the entanglement between the system s 
and an ancilla qubit a after the first and the second collision. This entanglement decreases when the 
system becomes more correlated with the reservoir and vice-versa, therefore, it identifies properly any 
backflow of information from fhe latter fo fhe former. We quanlify enfanglemenf by measuring fhe 
concurrence C ll32]| befween sysfem and ancilla and in Fig. [^) we plof ifs difference afler one and 
fwo collisions, C(2) — C(l), as a function of s. The values of C(l) and C(2) are obfained from fhe 
tomographic reconsfrucfion of fhe fwo-qubif densify mafrices and are plofed in fhe insef of Fig. [^). 
Fig. [^) shows a fransifion from a negative fo a posifive difference af around s = 0.3. Nofe fhaf posifive 
difference (C(2) > C(l)) means fhaf system and ancilla are more enfangled after fwo collisions fhan 
afler one which, in our model, idenlifies slrong non-Markovianily. As s decreases, fhe evolufion of fhe 
sysfem will mimic a sysfem fhaf gels more correlaled fo fhe reservoir after fhe second collision (and 
Iherefore less enfangled wifh fhe ancilla) which defines fhe regime of weak non-Markovianily where 
A 2 i’s non-CP character can only be evidenced by full lomography of fhe map ilself. Also note fhaf 
fhe Iheorelical curve predicls a disconlinuily in fhe derivative of C(2) — C(l) and a plateau for a range 
of values of s. Bolh behaviors are easily explained if we look al fhe individual behaviors of C(l,2) 
plotted in fhe insef of Fig. [^). Bolh enfanglemenls suddenly die ||33l, C(2) faster fhan C(l), bul 
C(l) remains zero for any larger value of s while C(2) evenfually recovers due to fhe environmental 
correlations. The regions of weak and strong non-Markovianity (blue and red regions, respectively) are 
also presented in Fig. [^); however, here there is no clear sign of the transition between one region to 
the other, so, if just the divisibility of the maps is calculated, there is no essential difference between 
these two types of non-Markovianity. The theoretical curves are computed assuming imperfections 
in the preparation of the initial state (C(0) ~ 0.975) and imperfect operations of the LC devices. We 
modelled the imperfections in the operations as Xexp(p) = FXpX-l- (1 —F)/2YpY -|- (1 —F)/2ZpZ 
and ’Zexpip) = F'Ep'E + (1 — F)/2YpY-|- (1 — F) /2XpX, and we considered F = 0.97. 

Discussion 

We have realized experimentally the collisional model proposed in Ref. [26] to investigate the non- 
Markovian dynamics of an open quantum system. We showed how the evolution of the same photonic 
system can transit from strong to weak non-Markovian evolution by varying only one parameter. 
This effect is caused by simply modulating the probability of the photon to undergo a rotation on its 
polarization state. As a result, a particular kind of non-Markovianity which is normally not spotted in 
other experiments is observed here. All non-Markovianity is caused solely by simulating a correlated 
reservoir. Finally, the fact that both regimes are produced by the same underlying physical mechanism 
explicitly shows that there is nothing necessarily fundamental about strong non-Markovian evolutions. 

Besides its intrinsic relevance on the fundamental side, the weak non Markovianity experimentally 
demonstrated in this work allows to envisage future important applications regarding for instance 
quantum control techniques and resolution enhancement in quantum metrology If9l[l0l. 

Methods 

The positivity and completely positive character of our map can be identified by the dynamical matrix 
!}{. For the intermediate map A 21 and error models previously explained, the reconstruction of gives 
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First, notice that, considering that our operations are almost perfect {F = 0.97), ?i 2 < 0 and, as already 
explained, the map A 21 is not completely positive. However, we would like to identify for which values of 
£ the map becomes positive. A map is positive if the dynamical matrix is block positive, i.e. Xi + Xj > 0, 
where Xi are the eigenvalues of ll^ . It is possible to show that the only inequality that does not satisfy 
this condition is Xq + X 2 - For the case of perfect operations {F = 1), the map will be positive only when 
(1-4£) >0. 

The witness that we used to track the backflow of information was the concurrence. For our maps, it 
is given by 

C(l) = Max[0,l-4£], (3) 

C(2) = Max[0,l+4£((F(3F-2) + 3)£-2)]. (4) 

For perfect operations {F = 1), we can see that the condition established for positivity of the map 
((1-4£) > 0) will definitely imply C(2) —C(l) <0. A monotonic decay of the concurrence will identify 
then when the map is positive, and on the other hand, non-positivity ((1 — 4£) < 0) implies an increase 
in the concurrence. However, when the operations are not perfect, there is a range of values of £ for 
which C(2) = C(l), and this witness will fail to identify the exact point of transition from weak to 
strong non-Markovianity. Nevertheless, it will still be true that if C(2) < C(l) the map is positive (weak 
non-Markovian) and if C(2) > C(l) the map is not positive (strong non-Markovian). 
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